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Positive Definite Solutions of the Nonlinear Matrix Equation 

X + A^X'^A = I 



Bin Zhou*^ Guang-Bin Cai-f James Lam^ 



Abstract 



^ ■ This paper is concerned with the positive definite solutions to the matrix equation X + A^X ^ A = I 

O ■ where X is the unknown and ^ is a given complex matrix. By introducing and studying a matrix 

operator on complex matrices, it is shown that the existence of positive definite solutions of this class of 
nonlinear matrix equations is equivalent to the existence of positive definite solutions of the nonlinear 
matrix equation W + B^W~^B — I which has been extensively studied in the literature, where B is a 
real matrix and is uniquely determined by A. It is also shown that if the considered nonlinear matrix 
■ equation has a positive definite solution, then it has the maximal and minimal solutions. Bounds of the 

positive definite solutions are also established in terms of matrix A. Finally some sufficient conditions 
and necessary conditions for the existence of positive definite solutions of the equations are also proposed. 
^ ■ Keywords: Bound of solutions; Complex matrix; Nonlinear matrix equation; Positive definite solutions. 



1 Introduction 

^ . Various kinds of matrix equations have received much attention in the hterature (see, for example, [1], [7], [5], 

^ ; [H], [12], [S], [m, [H, [IS], [13, dZ], [13, [H], [33, [SB], [33, [10], and the references therein). Especially, 

QQ . the problem of finding fixed points of the nonhnear matrix equation X + A* X^^ A = Q where A and Q > 

are given and X is unknown, has been extensively studied in the last two decades. The interest of studying 
such problem mainly relies on its applications in many fields such as the analysis of ladder networks [2], 
dynamic programming [33], control theory [22], stochastic filtering [1] and statistics [32] (see [3], [18] and 
PsJ . [in] for detailed introduction). Some generalized forms of the nonlinear matrix equation X + A*X~^A = Q 

have received much attention in recent years (see, for example, [6], [17], [24], and [30 ). 

Among the existing publications in the literature, two kinds of results can be found. The first kind of 
■ results concentrate on providing analytical conditions on the existence of positive definite solutions and their 

H ! corresponding properties. For example, the shorted operator theory was applied in (3] to study the existence 

of a positive definite solution; necessary and sufficient conditions in terms of symmetric factorizations of 
some rational matrix-valued function were derived in |18j for the existence of positive definite solutions; 
and some necessary and sufficient conditions were also derived [35j for the same problem in terms of some 
factorizations of the coefficient matrices. The other kind of results are mainly concerned with the numerical 
solutions of this class of nonlinear matrix equations. Basically, this can be accomplished via iterations 
including inversion- involved iterations [IH] , [H] and inversion- free iterations [2D] , [23] , [31] , [3S] • 

In the present paper, we consider a variation of this well-studied nonlinear equation. We study the nonlinear 
matrix equation X + A*X A = Q, which, as we will show in this paper, has totally different solutions 
from the solutions of X -I- A*X^^A — Q. In particular, we are interested in the existence of positive definite 
solutions of such kind of nonlinear matrix equations. Via some specific representations of complex matrices, 
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we are able to transform the equation X + A*X ^ A = Q into the equation W + B*W^^ B = P, where B and 
P are determined by A and Q, respectively. This allows us to study the original nonlinear matrix equations 
with the help of the existing results on the equation X + A*X^^A = Q. Other topics of this paper include 
the estimate of the bounds on the solutions, sufficient conditions, and necessary conditions on guaranteeing 
a positive definite solution. 

The rest of this paper is organized as follows. The problem formulation and some preliminary results to be 
used are given in Section [51 In Section [31 we present necessary and sufficient conditions for the existence of a 
positive definite solution of the considered nonlinear matrix equations. Both upper bounds and lower bounds 
of the solutions will be established in Section [4l while the necessary conditions and sufficient conditions 
guaranteeing a positive definite solution are given in Section [5l We will draw the conclusions of this paper 
in Section [6l 

Notation: In this paper, for a matrix A, we use A"^ ,A* , A, A (A) , det (A) , and p (A) to denote respec- 
tively the transpose, the conjugated transpose, the conjugate, the spectrum, the determinant, the 2-norm, 
and the spectral radius of A. Moreover, uj {A) = max{|z| : z — x*Ax, \\x\\ = 1} is the numerical radius of A. 
Finally, the symbol P > means that P is positive definite, /„ denotes an n x n identity matrix, denotes 
a zero matrix with appropriate dimensions, and j = \/— T. 

2 Problem Formulation and Preliminary Results 

We consider the following nonlinear matrix equation 

X + A*X'^A = Q (1) 

where Q E C"^" is a given positive definite matrix, A e C"^" is a given complex matrix, and X e C"^" is 
the unknown. In this paper, we are interested in the existence of positive definite solutions of this class of 
nonlinear matrix equations. 

Remark 1 Similar to 137} , we can also consider positive definite solutions of matrix equation 

X + A*X-^A = Q. (2) 

However, the positive definite solutions of (0) coincide with the positive definite solutions of equation (QJ) 
since X~^ =X~* =X~'\ 

Via a simple manipulation we can show the following result. 

Lemma 1 Let Q be a positive definite matrix. Then X is a solution of (Op if and only if Y = Q^^XQ^^ 
is a solution of the following nonlinear matrix equation 

I^=Y + A*qY''Aq, Aq^Q'^AQ-^. 

Therefore, without loss of generality, we assume hereafter that Q = /„ in ([T]). We point out that matrix 
equation (|T]) has solutions that are totally different from solutions of the following nonlinear matrix equation 

X + A*X-'^A^ In. (3) 

See the following example for illustration. 



Example 1 Consider a nonlinear matrix equation in the form of (Qp with Q = I2 and 



Then according to the results we will give later, we find the maximal positive definite solution of this equation 
as 



X+ = 



, -i - + 

8 ' SJ 2^8 



2 ' 8 ^ " 8 8J 

1,1. 1,1/6 



However, according to the results in 1 18^ . the maximal positive definite solution of equation (0) can be 
computed as 

It is clearly that Xj^ ^ X^. 



For a complex matrix A = A\+ A-i} G C"x™ where A\,A2 € R"^™, we denote the operators (•) and (•) 



as 



A^ = 



Ax -A2 
A2 Ai 



, A0 = 



A2 Ax 
Ai -A2 



It follows that both AP and A'^ are real matrices. For further use, we define two unitary matrices En and 



/„ 

/„ 



Pn. 



V2 



(4) 



Some basic properties of these matrices are collected as the following lemma whose proof is provided in 
Appendix B. 



Lemma 2 Let A e C"^™ and B e C™^'' he two given complex matrices. 

1. The following equalities are true. 

= (A)'' ^EnA^E^n, A^-i?„A^. 

2. Let Pn he defined in Then 



A"^ = Pn 



A 0_ 
A 



Consequently, A > (>)0 if and only if AP > (>)0. 

3. A £ C"'*" is a normal (unitary) matrix if and only if A^ is a real normal (unitary) matrix. 

4. For any A € C"^", there holds p (A) = p {A"^) and p {A^) = p^ {Ai) . 



5. The real matrix A"^ is normal if and only if A is con-normal, namely. A* A — AA* . 

6. For any A G C"^™, there holds \\A^\\ = = ||A|| 



7. IfA> 0, then {A'^Y ^ (a^) 



(5) 



Two matrices A,Be Qn-xn ^^^.g gg^j^ con-similar if there exists a nonsingular matrix S E C"X" such 

that S^^BS = A. The following lemma is borrowed from j4]. 

Lemma 3 Let Jk (A) denote a k x k .Jordan matrix whose diagonal elements are A. Then any matrix A € 
Qnx7i con-similar to a direct sum of hlocks of the form Jk (A) where X> or 

Ik 
Jk (A) 

where X < or 
mathrmlm (A) 7^ 0. 
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The set of complex numbers A appearing in the Jk (A) in the blocks of the the canonical form defined in 
Lemma [3] of a matrix A will be called the con-spectrum of A and denoted by co\{A) [3]. Moreover, the 
con-spectrum-radius of A will be denoted by 

cop {A) = max{|A| : A € coA {A)}. 

A property of the con-spectrum-radius is given in the following lemma whose proof will be presented in 
Appendix C. 

Lemma 4 Let A e C"'*" he a given matrix. Then cop {A) < (>) p {AX) < (>) 1. 
At the end of this section, we recall the well-known Schur complement. 
Lemma 5 Let matrix <1> be defined as 

^ ^ r <Z'ii <Pi2 ' 

[ ^1*2 ^22 

Then the following three statements are equivalent: 

1. <P>0. 

2. ^11 > and ^22 -~ ^i2^iY^i2 > 0. 

3. $22 > and - ^12^22^ ^12 > 0. 

3 Necessary and Sufficient Conditions 

In this section, we study necessary and sufficient conditions for the existence of a positive definite solution 
of the nonlinear matrix equation ([l}. Firstly we introduce an useful lemma. 

Lemma 6 Assume that A is nonsingular. Then X solves the nonlinear matrix equation (J^ if and only if 
Y ~ In — X solves the following nonlinear matrix equation 

In^Y + AY^^A*. (6) 

Proof. Let X be a solution of equation ([IJ, then A*X ^ A — In — X from which we get X ^ = 
A~* {In — X) A^^. Taking inverses on both sides gives X — A {In — X) A* which is equivalent to equation 
([6]) by setting In ~ X = Y. The converse can be shown similarly. ■ 

Our main result regarding the existence of positive definite solution of equation ([1]) is presented as follows. 

Theorem 1 The nonlinear matrix equation in |ip has a solution X > if and only if the following nonlinear 
matrix equation 

l2n^W+{A0fw-'A'>, (7) 
has a solution > 0. Moreover, the following two statements hold true: 

1. If the nonlinear matrix equation in flj) has a solution A" > 0, then it must have a maximal positive 
definite solution A+. Particularly, ifW+ denotes the maximal solution of the nonlinear matrix equation 
in then W+^X^, or 

+ 2 



In 



w+ 



iln 
In 



(8) 
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2. If A is nonsingular and the nonlinear matrix equation |ip has positive definite solution X > 0, then it 
must have a minimal positive definite solution X^. Particularly, if W- denotes the minimal solution 
of the nonlinear matrix equation in then W- = X'^, or 



(9) 



Proof. Let X > be a solution of equation ([T]). Taking (•) on both sides of equation ((TJ and using 

Lemma [5] gives 



A*X A 



= X"" + {A'^f ((Xf 



-1 



= X'' + {AOf {X'^y' AO, 
which indicates that W = X'^ >Oisa solution of equation 



(10) 



"-4=" Let equation ([7]) have a solution W > 0. Then it must have a maximal solution according to Lemma 
[7] in Appendix A. We denote such maximal solution by W+. In the following we show that there must exist 
a matrix Y > such that W+ — . According to Lemma [8] in Appendix A, we know that 



converges monotonically to W+, namely, < < M^fc+i < Wu < Iim^ > and 

lim Wk = W+> 0. 



(11) 
(12) 



(13) 



We show that, for any integer k > 0, there exists a matrix Yk > such that 

Wk = if, Vfc > 0. 

We show this by induction. Clearly, equation (IT3|) holds true for fc = by setting Yq — In- Assume that (fT3|) 
is true with k — s, say, there exists a Fs > such that Wg — Y^ . Then, for fc = s + 1, by applying Lemma 
m we have 



— I2n 



{A^fiYrr'A^ 



{A'^f (E„YyE. 



n^.^Ejy'A^ 



I,^-iA*f(Y 



'2ri 



A*Y, ^A 



where K^+i ^ I - A*Ys A. As Ws+i > 0, we know that Ys+i > also. Therefore, p3| is proved by 
induction. 

Hence, it follows from (IT2]) and (IT3l) that there exists a matrix Yqc > such that 



+ - lim Wk 

k^oc 



lim Yi, 



namely, 

which, by using a similar technique used in deriving ([TU)) . is equivalent to 
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Hence the nonlinear matrix equation ([1} has a solution X = Yao- 

Proof of Item 1: Since for any positive definite solution X of ([1]), there is a real positive definite solution 
of (O, we must have X^ < W+ — which indicates that X < Fqo- However, X = Y^o is also a solution 
of the nonlinear matrix equation Hence, the nonlinear matrix equation ([T]) has the maximal solution 
= Yoo ■ The relation W+ = then follows directly. 

Proof of Item 2: If the nonlinear matrix equation ^ has a positive definite solution X > Q, then, as X < /„, 
by Lemma ini the nonlinear matrix equation (|6|) also has a positive definite solution In — X, which, according 
to item 1 of this theorem, indicates that equation ([6]) must have a maximal positive definite solution Kf . 
Hence, by applying Lemma [S] again, X^ = /„ — is also a solution of equation ((T)). In fact, X^ is the 
minimal positive definite solution of equation ([ij. Otherwise assume that X^^ < X^ is a positive definite 
solution of (H]). Then 

= /„ — x^ > In — x_ = y+, 

is a positive definite solution of This is impossible since Y+ is the maximal positive definite solution of 
equation ([6|) by assumption. 

According to item 1 of this theorem, the maximal positive definite solution Y^ to equation ^ is related 
with = Y^ where is the maximal positive definite solution of 

l2n^Z + A0z-'{A0f . (14) 

As is nonsingular, by applying Lemma [6] again, the maximal positive definite solution Z4. to equation 
is related with = /2„ — Y_ where Y_ is the minimal positive definite solution of equation ([7]). Hence 
the minimal positive definite solution X_ to equation ([1]) satisfies 



= {In - Y+f = hn -Z+=Y^ 



The proof is finished. 



Remark 2 The nonlinear matrix equation in ^ is in the standard from of h39\) with Q = In (see Ap- 
pendix) which has been extensively studied in the literature. Therefore, by adopting the existing results on the 
nonlinear matrix equation we can get corresponding results on the original nonlinear matrix equation 

m- 

According to the proof of Theorem [l] we have the following result regarding iteration based numerical 
solution of the nonlinear matrix equation ([l}. 

Corollary 1 Assume that the nonlinear matrix equation has a positive definite solution. Denote the 
largest solution by X^. Then the iteration 



Wi 



k+l 



(15) 



converges and such that limfe_^oo Wk — A"^. Moreover, if ||A_|_"'^^|| < 1, then the iteration in HIS]) converges 
to A^ with at least a linear convergence rate, namely, there exists a k* > and a number < /i < 1 such 
that 



Wk+i - X) 



Wk-x': 



Vfc > k*. 



(16) 



Proof. The convergence of the iteration ([T5|) follows from the proof of Theorem [TJ So we need only to show 
(I16p . We use the idea found in [36 to prove the result. According to Theorem [1] we have X'^ = W+, the 
maximal solution to equation ([T]). Notice that 



X) 



\Wk+i-W+\\ 

hn {AOf W-'aO - (/2„ - (^0)^ W^'aO 

{A^f{W-'-W,-')A^\ 

{W+'A'>f [Wk - W+) W^^A'^ 
< ||W+M^|| ||T4^-M^|| \\Wk - W+\ 



(17) 



6 



-25 L 



10 



Figure 1: Numerical solution of the nonlinear matrix equation ([IJ via iteration ([TF 



On the other hand, by using Lemma [21 we can compute 



A"- 









En (^+^ 








1 N 









\X7^A\ 



Let ^ e (0, 1) be such that \\W^ ^^^\\ = ll^'^+^^ll < a/M < 1- As limfc_s.oo Wk ~ X'^, there exists a k* such 
that IIVFjT^A^II < y/Ji < l,Vfc > k*. Hence, the inequality in ([T7| reduces to (|T6|) immediately. The proof is 
finished. ■ 

We emphasize that the condition < 1 is only sufficient for guaranteeing the linear convergence of 

the iteration in (jlSp which converges as long as the nonlinear matrix equation ^ has a positive definite 
solution. By combining Lemma [6] and Corollary [1] we can also present a result regarding obtaining the 
minimal solution to the nonlinear matrix equation ([iJ. The details are omitted for brevity. 

Example 2 Consider the nonlinear matrix equation in Example{J\ By computation we have \\X^^A\\ ~ 
0.614 < 1. Then by Corollary [Jl we conclude that the corresponding iteration in il5\} converges to X"^ at 
least linearly. For illustration, the history of the iteration is recorded in Figure [7J From this figure we see 
that the convergence of the corresponding iteration in Iil5\) is indeed linear. Hence, the estimation in il6]) 
may be nonconservative. 

In the particular case that A is real, we can show the following result. 



Corollary 2 Suppose that A is real. If the nonlinear matrix equation LS9\} has a positive definite solution, 
then is real. Furthermore, if A is nonsingular, then X^ is also real. Hence in this case, equation (QJ) 
and the following nonlinear matrix equation 



In = X + A'X-'A, 
has the same maximal and minimal positive definite solutions. 
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Proof. Under the assumption of this corollary, W+ = X'^ and = X'2 are respectively the maximal and 
minimal solution of equation ([7]). Therefore, we need only to show that W+ and W- are in the form of 



Wi+ 
Wi+ 



Wi- 
Wi- 



Since W+ is the limit of the iteration in ([TT|) , we only need to show that 



Wik 
Wik 



Vfc > 0. 



(18) 



We show this by induction. Clearly, (1181) holds true for k ^ 0. Assume that it is true with k — s. Then, for 
A: = s + 1 , we can compute 



W. 



s+l 



In - AjWuAi 

/„ - ^^1^1,^1 







Therefore, (|18p is proved by induction. The case W- can be proved similarly. ■ 
Combining Theorem [T] and Lemma [10] gives the following corollary. 

Corollary 3 Suppose that A is invertible. Then the nonlinear matrix equation ^^ has a positive definite 
solution if and only if ui (^^) < ^■ 

Our next theorem presents some properties of the maximal and minimal positive definite solutions of the 
nonlinear matrix equation ([1}. 



Theorem 2 At 



that the nonlinear matrix equation in (QJ) has a solution X > 0. 



1. Let the maximal positive definite solution he X^. Then X^ is the unique solution such that X^ + 
sAX^ A is invertible for all s Cz {z : \z\ < 1} , or equivalently, 

cop [li+^A] < 1. (19) 



2. Assume further that A is nonsingular. Let the minimal positive definite solution be X^. Then X^ is 
the unique solution such that X- + sA X _ A* is invertible for all s G {z : \z\ > 1}, or equivalently, 



cop{X_ A*] >1 



(20) 



Proof. Proof of Item 1 : By Theorem [1] W+ — X'^ is the maximal solution of the nonlinear matrix equation 
([7|). Then according to Lemma [7] in appendix, Wj^ is the unique positive definite solution of equation ([7]) 
such that W+ + AA^ is nonsingular for all A € {z : |z| < 1}. Since P*£'„P„ = iJ„, we can compute 



W+ + AA<> =XX + XEnA"^ 



from which it follows that 



Pn 
Pn 
Pn 
Pn 



X+ 

x+ 

x+ _0 
Xh 

X+ 
Xh 

X+ XA 
XA X+ 



P* + XEnPn 

F xp:e„p, 

^XEn 

p* 



A 0_ 
A 

A 

A 0_ 
'A 



p^ 



det {W+ + XA^) = det {X+) det (^X+ - X^AX^^A 



(21) 
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Hence, det (W+ + AA^) is nonzero for all A G {z : |z| < 1} if and only if dct{X^ — \^AX_^_^A) is nonzero for 
all A e {z : |z| < 1}. Notice that 

VA e {z : \z\ < 1} ^ Vs = -A^ e {z : \z\ < 1} . 

The first conclusion then follows directly. Moreover, it follows from (|2T|) that 

det {W+ + XA'^) = det {X+) det (^/„ + sAX^^ AX^^'j det {X+) 
= (det iX+)f det (/„ + sAX;^AX+ , 

where s G {z : |z| < 1} . Hence det (W+ + XA^^ is nonzero for all A G {z : |z| < 1} if and only if the matrix 
AX AX 7^ has no poles ae{z:|z|>l} which is equivalent to 



1 > p (AXl'AX-') = p ( (x-'a) ( X-'A 



XI' a) (x-'a) 



By applying Lemma IH the above inequality is equivalent to (|19p . Assume that there exists another positive 

definite solution X^ such that cop (^X_^_' Aj < 1. As the above process is invertible, we can show that 

= Xf is such that W+ + XAO is nonsingular for allAG{z:|z|<l} which contradicts with Lemma [7] 
in appendix. 

Proof of Item 2: By Theorem [U Vr„ — X'^ is the minimal solution of the nonlinear matrix equation ([7|). 
Then according to Lemma [7] in appendix, W- is the unique positive definite solution of equation ([7]) such 
that W-+X{A0) is nonsingular for all A G {z : |z| > 1}. Similar to the proof of item 1 of this theorem, 
via some computation, we have 



W-+X{A^f ^ Pn 



X_ XA 
XA* X- 



p* 



The remaining is similar to the proof of item 1 and is omitted for brevity. The proof is complete. ■ 

4 Bound of the Positive Definite Solutions 

In this section, we study the bound of the positive definite solutions of the nonlinear matrix equation ([l}. 
Let {Hk}^^i be generated recursively as follows: 

m [Ml 

k is odd. 



Hi — In, Hk+1 — < 




(22) 



k is even. 



Here, if fc = 2, the zero matrices in H2 are obviously of zero dimension. Then we can prove the following 
result. 



Theorem 3 // the nonlinear matrix equation ([Ip has a positive definite solution, then > 0,Vfc > 1, and 
for any integer k > 1, there holds 



X > 



X > 






* 














A* 





A* 



= Sk, k is even, 



=^ Sk, k is odd, 



(23) 
(24) 
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in which, if k — 1, the zero matrices involved are of zero dimensions. Moreover, Sk is nondecreasing, namely, 

Sk+i>Sk, Vfc>l. (25) 



Proof. Since X — In — A*X A > 0, via a Schur implement, we have 

n ^ 4! 1 - [ ^" ^* 
[A X \ [A In- A*X- 

which indicates that H2 > 0. Applying another Schur complement on the above inequality gives 



















. A* 







< 


H2 






. A* 






H2 











[[0 A] 


X 




_[0 A] 


In- 




-'a 





A* 



X 





A* 



which indicates that > 0. By using Schur complement again, the above inequality can be equivalently 
rewritten as 



< 


H3 






A* 






Hi 






A* 






[[0 A] 


X 




[[0 A] 


In- 


-Tx 


-'A 













X-' 




. 3* 





which indicates that II4 > 0. Repeating the above process gives Hk > 0, Vfc > 1 
On the other hand, from the above development, we see that, if k is even, then 

> 0, 



and if k if odd, then 




>0, 



where, if fc = 1, the zero matrices are of zero dimension. We notice that the above two inequalities are 
respectively equivalent to and (|24l) via Schur complements. 



We finally show (pS)) . We only prove the case that k is odd since the case that k is even can be proven 
similarly. As A: + 1 is even, we have 



Sk+i — 





T 


A* 



k+l 





A* 



[0 A] 





A* 



In 





A* 



(26) 



We notice that, if k is odd, then 
Hk 



[0 A] 





A* 



A 


Vii 



[A* ] 



Hu 







■ A ' 


(Hk 








V12 

[A* ] 



10 



where Vn and V12 are two matrices that are of no concern. With this, the relation in (PB)) can be continued 



as 



Sk+i 





A* 



A 




[A* ] 



> 





A* 


T 



HI' 



A 




[A* ] 





A* 


A* 





A* 



= Sk. 

In the above development, the zero matrices involved are of zero dimension if fc = 1. The proof is finished. 
■ 

Since Sk is bounded above, there must exist a positive definite matrix 5*00 < In such that 

hm Sk = Soo- 

Obviously, we have X > Soo- It is not clear whether 5*00 — X^. Moreover, it is clear that the larger the fc, 
the better the Sk gives a lower bound of X. However, large k may lead to numerical problems. By choosing 
some special values in k, the following corollary can be obtained. 



Corollary 4 // the nonlinear matrix equation flj) has a positive definite solution, then A satisfies 



In > A A* + A* A. 

Moreover, the solution X satisfies the following inequalities 



X > Si = AA*, 

X > S2=A{In~AA* 
X > S3^A(ln-A {In - AA* 



A' 



(27) 

(28) 
(29) 

(30) 



Proof. Let fc = 2. We get from Theorem |3] that H2 > 0, which, via a Schur complement, is equivalent to 
/„ > AA* . Similarly, by letting fc = 3, we get from > that 



which, via a Schur complement, implies 

In-[0 A] 





" /„ 


A* 





< 


A 


In 


A* 







A 


In 



'In A* ' 


-1 





A /„ 







= I„-A{In-AA*)-^A^>0. 



By applying Schur complement again, the above inequality is equivalent to 

/„ A 
A^ In - AA* 



> 0. 



However, another Schur complement indicates that the above inequality implies A*A + AA* < In 

Now if we set fc = 1, we get from ^ that X > A A* = Si (or X > AA'^ = 5*1) which is 
k — 2, we get from ((23)) that 



if we set 



X > S2^[0 A] 
= [0 A] 



A{In-AA*)-'A'^, 







A* ' 


-1 







A 


In . 






' Vii 




V12 










. ^2 


{In 


-AA*y 


1 









A* 
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which is (P^ . Here V^j denotes the elements that are of no concern. Moreover, if we set fc = 3, we know 
from ^ that 



X> S3 





A* 



= [0 A] 



A 







.A* 




In A* 





-1 


A In 


A* 




A 


In 




Vii 






Vl2 ( 


I 71 A (/„ 





A* 
V12 





T 



A(^In-A{ln-AA*) 'a*^ 'a^, 



which is just (|30p . The proof is finished. ■ 

With the help of Lemma [6] and Theorem [3l we can also obtain upper bounds of X. To this end, we let 
{^felfcli be generated as where ^* is replaced with A, namely. 



Gk 



Gl — In, Gk+l — < 



[OA*] 
Gk 





A 



A 



In 



A 



k is odd 



fc is even. 



Again, if fc = 2, the zero matrices in G2 are obviously of zero dimension. 

Theorem 4 Assume that A is nonsingular and the nonlinear matrix equation in (QJ) has a positive definite 
solution. Then Gk > 0,Vfc > 1, and for any integer fc > 1, there holds 



and 



X < In 



X<In- 





A 





A 



Gl. 



Gk' 





A 





A 



= Rk, fc is even, 



= Rk, fc is odd. 



(31) 



(32) 



in which, ifk — 1, the zero matrices involved are of zero dimension. Moreover, Rk is non-increasing, namely, 

Rk+i < Rk, Vfc > 1. (33) 



Proof. Since A is nonsingular, by Lemma IHl equation ([1]) is equivalent to equation ^ with Y = I — X. 
Applying Theorem [3] on equation (|6]) gives Gk > 0, fc > 1 and 



— Sk, fc is even, 



— S'/., fc is odd. 



which are respectively equivalent to (pij) and ([5^ by using Y = I—X. Finally, Rk = /„ — S*^ is non- increasing 
because S*^ is nondecreasing according to Theorem [S] The proof is finished. ■ 

It is also clear that the larger the fc, the better the Rk gives an upper bound of X. However, large fc may 
lead to numerical problems. Choosing some special values in fc gives the following corollary. 



Y > 




_ A 


Gk' 




A 


Y > 




A 


* 

Gk 




A 
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Corollary 5 // the nonlinear matrix equation in |ip has a positive definite solution and A is nonsingular, 
then A satisfies 

In > A*A + AA*. 
Moreover, the solution X satisfies the following inequalities 

X < i?i = /„ - A* A, 



X < i?2 = /„ - A* {In - A* A) A, 

X < i?3 - In - A* (in - A^ {In - A* A)'^ A) A 



5 Sufficient Conditions and Necessary Conditions 

In this section, we present some necessary conditions and sufficient conditions for the existence of a positive 
definite solutions of the nonlinear matrix equation ([l}. 



Theorem 5 // the nonlinear matrix equation ^j) has a positive definite solution, then p {AA^ < j, \\A\\ < 1 

and 

p[[a±A*)[a±T))<1, (34) 

which can be equivalently rewritten as cop [Azt A ) < 1. 



Proof. It follows from Theorem [T] and Lemma [TT] in appendix that if the nonlinear matrix equation ([l} has 
a positive definite solution, then p {A'^) < i, ||^^|| < 1 and 



p(ao±(a^)^) <1. 



(35) 



Clearly, by applying Lemma [21 p (A^) < ^ is equivalent to p {AA) < j and ||A^|| < 1 is equivalent to 
\\A\\ < 1 (we point out that ||A|| < 1 also follows directly from (P7)) ). We next show that ([55)1 is equivalent 
to p4)) . By virtue of Lemma [2] and in view of P*i?„P„ — En, we get 



[AO±{A^f)=p[EnA'^±{A'^fEr 



= P [EnPn 
= p {p*EnP, 

^p(e, 



A 0_ 
A 



A 0_ 
A 



^n ^ 
± 



A 0_ 
A 



± 



A* 
'A 

A* 
A 
A* 
A 



P* E 



P*E P 



Er. 



-, 2N 



A±A* ' 
A±A* 

A±A* 
A±T 

{A±A*) (yl±T 


= p^[[a±a*) {a±a*)), 

which is the desired relation. Finally, the equivalence between (p4|) and cop (^A i ^ ) 
Lemma 21 The proof is finished. ■ 



A±A*) {A±A*) 



< 1 follows from 
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Theorem 6 The nonlinear matrix equation ([7P has a positive definite solution provided 

U\\ < \. (36) 

Moreover, if A is con-normal, then the nonlinear matrix equation (QP has a positive definite solution if and 
only if A satisfies fiSb]] . In this case, the maximal solution is given by 

^+ = i (/« + (/„ -4^M)^) . (37) 

// A is further assumed to be nonsingular, then the minimal solution can be expressed as 

X- (/«-(/« -4A*A)^). (38) 

Proof. If ||A|| < i, then by Lemma[21 we know that ||A^|| < \, which, by usmg Lemma WI[ indicates that 
equation ([7]) has positive definite solution. This is further equivalent to the existence of positive definite 
solution of equation ([T]) in view of Theorem [T] The case that A is con- normal can be shown similarly. We 
next show ((37)) . Notice that, according to Lemma fT2l the maximal solution of equation ([7]) is given as 

= \{l2n+[hn-^A'^f ElE^A'^Y^ 

^\{l.n+{hn-^{A^fA^Y^ 

= \[l2n+(l2n-A{A*Afy^ 

= -(/„ + (/„-4AM)^j , 

which, according to Theorem [Tl implies (j37p . The equation (1381) can be shown similarly. The proof is done. 
■ 

6 Conclusion 

This paper has studied the existence of a positive definite solutions to the nonlinear matrix equation X + 
A* X A — Q. With the help of some operators associated with complex matrices, we have shown that 
the existence of a positive definite solution of this type of nonlinear matrix equation is equivalent to the 
existence of a positive definite solution of a nonlinear matrix equation in the form of Vl^ + B"^W~^B — I 
where B is real and is determined by A. Since the later nonlinear matrix equation has been well studied in 
the literature, properties of the original nonlinear matrix equations can be established based on the existing 
results on the transformed nonlinear matrix equations. Moreover, with the help of Schur complement, we 
have shown in this paper some upper bounds and lower bounds on the solutions to the nonlinear matrix 
equations. Simultaneously, some easily tested sufficient conditions and necessary conditions for the existence 
of positive definite solution of the nonlinear equations have also been established. We point out that, by 
combining the results obtained in this paper and the existing results on numerical computation of solutions 
to the standard nonlinear matrix equation W + B"^W~^B = I, numerical reliable algorithms can be built 
for computing the positive definite solutions to the original nonlinear matrix equation, which is currently 
under study. 



14 



Appendix 



A: Solutions of Matrix Equation X + A*X-^A = Q 

In this subsection, we recall some basic results regarding positive definite solutions of the following matrix 
equation 

X + A*X-^A = Q. (39) 

Lemma 7 (Theorem 3.4 in J18^ ) Suppose that Q > and assume that he nonlinear matrix equation 139\) 
has a positive definite solution. Then it has a maximal and minimal solution and respectively. 
Moreover, X^ is the unique solution for which X + XA is invertible for all \X\ < 1, while X- is the unique 
solution for which X + XA* is invertible for all \X\ > 1. 

Lemma 8 (Algorithm 4.1 in JT^I) Suppose that Q = I. If the nonlinear matrix equation \39(l has a positive 
definite solution, then the iteration 

Xk+l = In — A*X^, "^A, Xq = /„, 
converges to the maximal solution X+ , namely, limfe^oo Xk = X+ . 

Lemma 9 (Theorem 8.1 in JlSf) Suppose that Q — I and A is real. If the nonlinear matrix equation I139\) 
has a positive definite solution, then Xj^ is real. Furthermore, if A is nonsingular, then X^ is also real. 

Lemma 10 (Theorem 5.1 in \lSf ) Suppose that A is invertible. Then the nonlinear matrix equation i39\) 
has a positive definite solution if and only if u) (A) < i . 

Lemma 11 (Theorem 7 in 119] and Theorem 3.1 in f35}j) If the nonlinear matrix equation i39\) has a positive 
definite solution, then p{A) < ^,\\A\\ < 1 and p {A ± A*) < 1. 

Lemma 12 (Theorem 11 and Theorem 13 in '19]) The nonlinear matrix equation Ii39\) has a positive definite 
solution provided \\A\\ < i. Moreover, if A is normal, then the nonlinear matrix equation i39\) has a positive 
definite solution if and only if \\A\\ < \. In this case, the maximal solution is given by 

^+-i(/„ + (/n-4ATA)^). 

Furthermore, if A is nonsingular, then 

X- = i(/„-(/n-4ATA)^). 

B: Proof of Lemma [2] 

Proof of Item 1: These equalities can be verified directly by definition. 
Proof of Item 2. This result follows from Lemma 9 in jS?] . 

Proof of Item 3: We need only to show that A is a normal matrix if and only if A^ is a real normal matrix 
since unitary matrix is a special case of normal matrix. If A^ is a real normal matrix, then (^'^) = 

A'' [A"^)^ . However, A"^ {A"^)^ = {AA*f and {A"^)^ A"^ ^ {A*Af . Hence we have {AA*f ^ {A*Af 
and consequently AA* = A* A, that is, A is a normal matrix. The converse can be shown similarly. 
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Proof of Item J^: From item 2 of this lemma, we obtain 

P (A^) - P (p. 



A 0_ 
A 



' A 


■ 







A 


) 



max{p(A),p(A)} 
p{A). 



Similarly, we can compute 



A 0_ 
A 



p (AO) = p (i?„A^) = p ( E,,P, 
= p (^P*E^P, 

^p{e, 

= P 



A 0_ 
A 



P* P 



P' 



A 0_ 
A 

A 
A 



AA 



^ 
A 



A 
A 



= P^ (AA) 



Proof of Item 5: By definition, we can compute 

= {A'^f A"" = iA*f A"" ^ iA*Af , 



and similarly, 



AO {A^f = Er.A'^ {EnA^^f = (A^)^ 



Clearly, is a normal matrix if and only if (y4^)'^A^ = A^ (A^)"^, which is equivalent to (AA*)^ = 
(AM)"^, namely, A*A^ XF. 

Proof of Item 6: By using item 1 of this lemma, we obtain 



U^ll = = lU^II = 



Pn 



A 0_ 

A 



p* 



A 0_ 

A 



= WM 



Proof of Item 7: Let U he a. unitary matrix such that A = UDU* where D is a real diagonal positive 
semi-definite matrix. Then 



Q\ 2 



[UDW 



(7^ 



D 
£> 



Di 
£15 



{u''f=(uD'^U 



C: Proof of Lemma [4] 

We only show the case "<". According to the results in [4], we know that 
(1). for A > 0, A e coA (A) ^ A^ g A (AA) . 
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(2). for A < 0, A e coA (A) ^ X G X (AA) 



(3). for Im (A) ^ 0, A e coA {A)^X,XgX (AA) . 

Let s be an arbitrary eigenvalue of AA. Then |s| < 1. Consider three cases. Case 1: 1 > s > 0. Then it 
follows that A = y/s G coA (A) and hence |A| < 1. Case 2: — 1 < s < 0. Then we see that A = s e coA (A) 
and hence |A| < 1. Case 3: Im (A) 0. In this case, we see that either A = s e coA {A) or A = s S coA {A) . 
In both cases, there holds |A| = \s\ < 1. The proof is completed. 
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